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The Atacama pathfinder experiment Sunyaev-Zel'dovich (APEX-SZ) instrument is a millimeter-wave cryo- 
genic receiver designed to observe galaxy clusters via the Sunyaev-Zel'dovich effect from the 12 m APEX 
telescope on the Atacama plateau in Chile. The receiver contains a focal plane of 280 superconducting 
transition-edge sensor (TES) bolometers instrumented with a frequency-domain multiplexed readout sys- 
tem. The bolometers are cooled to 280 mK via a three-stage helium sorption refrigerator and a mechanical 
pulse-tube cooler. Three warm mirrors, two 4 K lenses, and a horn array couple the TES bolometers to the 
telescope. APEX-SZ observes in a single frequency band at 150 GHz with 1' angular resolution and a 22' field- 
of-view, all well suited for cluster mapping. The APEX-SZ receiver has played a key role in the introduction 
of several new technologies including TES bolometers, the frequency-domain multiplexed readout, and the 
use of a pulse-tube cooler with bolometers. As a result of these new technologies, the instrument has a higher 
instantaneous sensitivity and covers a larger field-of-view than earlier generations of Sunyaev-Zel'dovich in- 
struments. The TES bolometers have a median sensitivity of 890 /xKcmb y/s (NEy of 3.5 x 10 -4 y/s). We have 
also demonstrated upgraded detectors with improved sensitivity of 530 iiKqmb\/s (NEy of 2.2xl0 -4 \/s). 
Since its commissioning in April 2007, APEX-SZ has been used to map 48 clusters. We describe the design 
of the receiver and its performance when installed on the APEX telescope. 
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I. INTRODUCTION 
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Galaxy clusters are the most massive gravitationally 
collapsed objects in the Universe, offering a unique probe 
into the composition and evolution of the universe (see, 
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e.g., Refs. [T]and [3]). Clusters of galaxies formed as dark 
energy became the dominant component of energy den- 
sity in the universe. The growth of structure from this 
period on is sensitive to the dark energy density and its 
equation of state. Therefore, measuring the evolution of 
cluster number density is a powerful method to constrain 
these parameters as well as the matter power spectrum 
normalization and the mass density of the universeP In 
addition, galaxy clusters are important laboratories for 
studying various astrophysical phenomena in great de- 
tail, like heat transport and instabilities in magnetized 
plasmas,^ evolution of galaxies in dense environments^ 
and cosmic circulation of heavy elements P 

The bulk of the baryonic mass in galaxy clusters re- 
sides in the form of a hot, low-density ionized gas filling 
the space between its constituent galaxies. This intra- 
clustcr plasma is observable at mm-wavelengths via the 
Sunyaev-Zel'dovich effect (SZE)Pa distortion of the cos- 
mic microwave background (CMB) blackbody spectrum 
arising from inverse Compton scattering of ~ 1% of the 
CMB photons by intracluster electrons. At 150 GHz, 
the SZE blackbody distortion is visible as an intensity 
decrement in the CMB. The surface brightness of the 
SZE signal is proportional to the integrated gas pressure 
through the cluster. The signal is independent of cluster 
redshift, since the SZE produces a fractional change in 
the CMB intensity which remains constant as the light 
is redshifted between cluster and observer. In contrast, 
x-ray and optical surface brightness dim with increasing 
redshift. The redshift independence makes the SZE a 
uniquely sensitive method for discovering and observing 
distant clusters (see, e.g., Refs. [Hand [9]). 

The Atacama pathfinder experiment Sunyaev- 
Zel'dovich (APEX-SZ) instrument is a millimeter-wave 
receiver for the APEX telescope designed to make 
sensitive measurements of the SZE in galaxy clusters.^ 
The APEX-SZ focal plane contains a 280 element 
superconducting transition-edge sensor (TES) bolometer 
array operating at 280 mK. The receiver-telescope com- 
bination has 1' resolution at its observation frequency 
of 150 GHz and a 22' field-of-view. Specifications and 
measured performance for APEX-SZ are summarized in 
Table [H 

The APEX-SZ detector array comprises six triangular 
sub-arrays with 55 bolometers each. The receiver was 
first installed on APEX with a single sub-array in De- 
cember 2005 for instrument characterization. We rede- 
ployed the receiver with a full focal plane array of 280 
multiplexed detectors in April 2007. Since then, we have 
accrued 875 h of observations in seven observing runs, 
using observing time allotted to German and Swedish 
institutions. Over this period, APEX-SZ has used two 
types of detectors with different absorption efficiencies 
and thermal links. We refer to these as type-1 and type-2; 
Sec. MI Cl details the differences between the two detector 
designs. 

APEX-SZ focuses on wide-field imaging of individual 
clusters to probe the relationship between SZE flux and 



Optics 

Primary mirror physical diameter (m) 
Primary mirror illumination diameter (m) 
Field-of-view (amin) 
Beam FWHM (asec) 



Detectors 
Bolometers 
Live channels 

Performance 
Band center (GHz) 
Bandwidth (Au) (GHz) 
Cumulative efficiency (ritot) 
NET (fiKcMB^) 
NEy (xl(T 4 v 7 *) 



12 

8 
22 
58 



280 
170-180 



Type-1 

151 
24.5 
0.31 

890 

3.5 



Type-2 

154 
33.5 
0.36 

530 

2.2 



TABLE I. APEX-SZ specifications and measured perfor- 
mance. The number of detectors typically used in data anal- 



ysis are listed as "live channels" (see Sec. VII A I. Bolometer 
performance values are listed for both type-1 and type-2 de- 
tectors (Sec. \Vll\ . The NET and NEy are listed for a single 
detector and refers to the median performance of the array 

(Sec.[ynH|. 



cluster mass at varying redshift. It is one of a number of 
millimeter-wave bolometer instruments which have come 
online recently The Large APEX Bolometer Camera 
(LABOCA; Ref. Q3J and AzTECP both observe from 
the Atacama at 345 and 270 GHz, respectively. Obser- 
vations of clusters at 150 GHz and higher frequencies 
are complementary, providing additional insight to clus- 
ter dynamics and constraining contamination from other 
sources. The South Pole Telescope (SPT; Ref. QJ and 
the Atacama Cosmology Telescope (ACT; Ref. IT4"|) are 
two dedicated SZE survey instruments which have dis- 
covered hundreds of clusters.^ APEX-SZ is more suited 
for targeted observation of clusters for multi-frequency 
study. The Atacama location enables SZE study of ex- 
isting x-ray data fields which is limited from the South 
Pole. The ACT is designed for surveying and does not 
have an elevation drive for pointing at individual clus- 
ters. The SZE-mass relationship measured by APEX-SZ 
will be critical for developing the cosmological constraints 
from the results of these surveys. 

SZE observations are functionally CMB anisotropy 
measurements at high sensitivity and resolution, and 
APEX-SZ represents a generational leap from the bolo- 
metric receivers that measured the primary anisotropics 
at degree scales. These receivers (such as the Mil- 
limeter wave Anisotropy experiment IMaging Array 
(MAXIMA; Ref. HU) and Balloon Observations Of 
Millimetric Extragalactic Radiation and Geophysics 
(BOOMERanG; RefP) each with 16 detectors) used 
neutron-transmutation-doped (NTD) bolometers. The 
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bolometers were housed in individual integrating cavities 
with separate feed horns and filters. Developed later, 
the BolocauP^ receiver employs a focal plane array of 
144 NTD bolometers, in which the absorbers are all built 
lithographically on a single silicon wafer. The focal plane 
includes a feed horn array and a single, large filter. 

In constructing APEX-SZ, we have expended much 
effort in developing robust, scalable technologies — 
notably TES bolometers and frequency-domain multi- 
plexed SQUID readout compatible with the pulse tube 
cooler — which represent a significant advance and en- 
able a large increase in sensitivity compared with previ- 
ous bolometric receivers. The bolometers are fabricated 
on silicon wafers with thin-film deposition and optical 
lithography. The focal plane and readout are assem- 
bled from several modular components enabling scaling 
to larger arrays. These components have been used in 
SPT, which started observations in 2007 with 960 de- 
tectors, and POLARBEARp^ an upcoming Berkeley CMB 
polarization experiment with over 1200 detectors. The 
analog multiplexed SQUID readout developed for APEX- 
SZ has evolved into a next-generation digital multiplexed 
readout system described in Ref. [201 The digital multi- 
plexer will be used in Polarbear, SPTpolpS and the 
E and B Experiment (EBEX; Ref. 122). 

This paper describes the design and performance of 
APEX-SZ. The structure of the paper is as follows: the 
APEX telescope and site are discussed in Sec. [TTJ optics 
in Sec. TTTJ TES detectors in Sec. |IV| multiplexed readout 
in Sec. VJ cryostat and cooling systems in Sec. |VI[ sy stem 
performance in Sec. VII observations in Sec. VIII[ and 
summary and conclusions in Sec. |IX| 



II. APEX 

A. Telescope 

The APEX telescope (shown in Fig. [IJ is a 12-meter 
diameter on-axis Cassegrain telescope based on the de- 
sign of the VERTEX Atacama Large Millimeter Array 
(ALMA) prototype, but modified to include two Nas- 
myth cabinsP^ the APEX was commissioned by the 
Max Planck Institut fur Radioastronomie, the European 
Southern Observatory, and the Swedish Onsala Space 
Observatory for use with bolometric and heterodyne re- 
ceivers and began operations in September 2005. 

The primary mirror consists of 264 panels aligned via 
radio holography to form a surface with 18 fim rms accu- 
racy allowing observations up to 1.5 THz. The secondary 
is 0.75 m in diameter and is mounted on a hexapod for 
optical adjustment. It focuses through a 0.75 m hole in 
the primary, which limits the field-of-view (FOV) to 22'. 
This FOV allows APEX-SZ to map the large angular size 
of low-redshift clusters. The primary and secondary mir- 
rors have been etched to scatter infrared and optical light, 
enabling daylight observations. The absolute pointing ac- 
curacy is measured to be 2" rms. The APEX secondary 




FIG. 1. The 12 m APEX telescope on Llano de Chajnantor 
at 5107 m altitude. 



support legs are made from carbon-fiber reinforced plas- 
tic and tapered to minimize scattering, and the large size 
of the primary mirror allows it to be under-illuminated 
while still achieving arcminute resolution. These two fea- 
tures are important for minimizing ground contamination 
since the telescope lacks a ground shield. 



B. Atacama Site 

The APEX site is at 5107 m altitude on Llano de Cha- 
jnantor in the Atacama Desert in northern Chile. The 
Atacama is one of the driest areas on the Earth and con- 
sidered, along with the South Pole and Mauna Kea, to 
be one of the premier locations for millimeter and sub- 
millimeter wave astronomy. 

The APEX is located near the ALMA site and within 
a few kilometers of the sites of a number of past (TOCO; 
Ref. [2j and CBI; Ref. [25]) and current (ACT; Ref. M 
and QUIET; Ref. |26| CMB experiments. The ALMA 
site testing campaign made atmospheric measurements 
at 225 GHz over several years and found atmospheric 
opacity at Chajnantor and the South Pole to be compa- 
rable with median zenith optical depths of t 2 25 = 0.061 
and t 2 25 = 0.053, respectively.^ While the opacities are 
comparable, the relative contribution from water vapor 
is greater in the Atacama desert than at the South Pole. 
Spatial variation of emission from water vapor results in 
higher fluctuation power at ChajnantorPSl Diurnal and 
annual variations in weather at Chajnantor can have 
a significant impact on the observed atmospheric noise. 
For much of the year, conditions arc excellent, but dur- 
ing much of the Austral summer, the "Bolivian winter" 
weather phenomenon brings 2-3 months of increased at- 
mospheric opacity and frequent precipitationPsl 

In addition to favorable observing conditions, the low 
latitude of the APEX site (23:00:20.8S, 67:45:33.0W) 
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allows observations of a large fraction of the celes- 
tial sphere. A large number of fields with rich multi- 
frequency data sets not accessible from the South Pole are 
observable from Chajnantor. Overlapping observations 
at several wavelengths enable a more complete study and 
characterization of clusters. For example, cluster ob- 
servations within the x-ray Multi-Mirror Mission-Large 
Scale Survey (XMM-LSS) field offer immediate compar- 
ison with x-ray data, and the Cosmological Evolution 
Survey (COSMOS) field provides a rich set of weak lens- 
ing, x-ray, and sub-millimeter data to cross-calibrate and 
characterize foregrounds for SZE observations. 



III. OPTICS 

A. Tertiary Optics 

The APEX-SZ tertiary optics couple the telescope 
beam to the focal plane, which is an array of feed-horn 
coupled bolometers. The optical system has diffraction 
limited performance over the full 22' FOV of APEX, and 
the focal plane is flat and telecentric, with //2.3. An 
additional constraint on the optical system is that the 
pulse-tube cooler cooling power varies with orientation 
(Sec. ~ 



VI B ) . The cryostat optical axis is tilted 30° from 



the telescope optical axis to keep the pulse-tube cooler 
within 30° of vertical over a 60° range in elevation angle. 



Cassegrain 
focus 



Folding flat 
mirror / 




Tertiary mirror 



Reimaging 
mirror 



FIG. 2. The APEX-SZ reimaging optics design. Each bundle 
of rays represents the beam from a single detector. 

Figure [2] illustrates the APEX-SZ tertiary optical sys- 
tem, which comprises three ambient temperature mir- 
rors, a 4 K Lyot stop, and two 4 K high-density polyethy- 
lene (HDPE) lenses. The tertiary mirror is a 1440 mm 



diameter off-axis paraboloid on the cabin floor. The ter- 
tiary is followed by a 620 mm flat (quaternary) and a 
980 mm x 1160 mm ellipsoidal mirror (quinary) that 
creates an image of the primary inside the cryostat. All 
three mirrors are pure conic sections. The cryostat win- 
dow aperture is 150 mm in diameter and is made from a 
50 mm thick layer of Zotefoam PPA30^ epoxied into 
an aluminum ring. Zotefoam is a nitrogen-expanded 
polypropylene foam, which has low scattering and ab- 
sorption at millimeter wavelengths. We measure the 
transmittance to be greater than 99% averaged across 
the instrument frequency response band. 

The 72 mm diameter 4 K Lyot stop is located at the 
primary image, and it truncates the beam at ~ 0.7 of the 
primary diameter to prevent spillover and ground con- 
tamination without significantly increasing the bolome- 
ter optical load. The edge taper at the Lyot stop is 
—4.5 dB, corresponding to a fractional spillover power 
of 0.34 for the D = 1.4/A diameter conical horns, where 
/ is the focal ratio and A is the optical wavelength. The 
V FWHM beams on the sky are the far-field diffraction 
pattern of this illumination profile. The horn diame- 
ter, and correspondingly, the Lyot stop edge taper and 
spillover power, were chosen to maximize array mapping 
speed given two design constraints: a maximum 22' FOV 
and a maximum of 280 readout channels. (An unlimited 
FOV favors larger horns, whereas an unlimited number of 
readout channels favors smaller horns.) The Lyot stop is 
blackened with a 7 mm thick layer of Eccosorb MF-117p^ 
which is machined with triangular grooves to reduce re- 
flections. 

All optical elements on the sky-side of the Lyot stop 
are sized wherever possible to keep fractional spillover 
power below 0.001 for a pixel on the edge of the array, 
see Table [TlJ This design minimizes both optical load- 
ing and truncation of diffractive sidelobes generated by 
the Lyot stop. Such truncation would degrade the qual- 
ity of the primary aperture stop. To assess the effect 
of finite optical element size on the beams and primary 
aperture stop, we modeled the beam propagation using 
multi-moded Gaussian beam propagation formalism and 
the ZEMAX physical optics propagation software.^ Both 
methods gave similar results, and indicated that our op- 
tical element sizing was adequate. 

Two 4 K high-density polyethylene (HDPE) lenses af- 
ter the Lyot provide the final beam shaping. To max- 
imize transmittance and minimize scattered light in- 
side the cryostat, the HDPE lenses have machined an- 
tireflection surfaces with closely spaced, circumferential 
grooves which have triangular cross-sections 0.56 mm 
wide and 0.64 mm deepP2 In effect, these grooves pro- 
duce a graded index surface coating between the vac- 
uum and the HDPE, which has an index of refraction 
n = 1.567 at 4 KP^ We calculate that the reflection 
loss at each lens surface should be less than 2% aver- 
aged across the band. The imaging lens axis is tilted 
with respect to the optical axis to compensate for the 
image tilt induced by the off-axis mirrors and create a 
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flat (telecentric) focal plane. 



B. Filters and Band 



APEX-SZ observes in a single band centered at 
150 GHz near the maximum intensity decrement of the 
SZE at ~ 130 GHz. The measured observation pass- 
band (shown in Fig. [3]) is selected to minimize contamina- 
tion from other astronomical signals and avoid the atmo- 
spheric molecular absorption lines for oxygen at 118 GHz 
and water vapor at 183 GHz. The pass-band is deter- 
mined by a series of low-pass (LP) optical filters^ (de- 
scribed below) , cylindrical waveguides and the bolometer 
absorption cavity (described in Sec. |III C 1. The broad 
water-vapor line at 183 GHz contributes roughly 3% 
to t he tot al background load on the bolometers (see 
SeclVHGI). 
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FIG. 3. Calculated atmospheric transmittance at APEX for 
2 mm (solid line) and 1 mm (thin solid line) PWV (precip- 
itable water vapor) at 55° elevation plotted with the measured 
APEX-SZ observation band (for type-2 detectors) (dashed 
line). The atmospheric transmittance data are from the 
APEX transmittance calculator based on the ATM modelP^ 

A series of LP optical filters serves to minimize radia- 
tive loading on the interior of the cryostat as well as de- 
fine the observation band. The LP filters consist of layers 
of capacitive, resonant metal-mesh in polyethylene. The 
filter scheme is shown in Fig. |4j The first element is a 
thin single layer LP filter, with metal-mesh on both sides 
of a polyethylene membrane, called an infrared shader. 
It is designed for a slow cut-off with minimal scattering 
loss. The ~12 THz (400 cm" 1 ) LP shader reduces the 
radiative load on subsequent filters. We minimize load- 
ing from 300 K emission on the 4 K stage with 3.0 THz 
(100 cm- 1 ) and 2.4 THz (80 cm" 1 ) multilayer LP filters 
at 60 K. Similarly, a 255 GHz (8.5 cm" 1 ) multilayer LP 
filter at the front of the 4 K radiation shell minimizes ra- 



diative loading from 60 K on the millikelvin stage. The 
series of filters with staggered cutoffs also serves to block 
harmonic leaks that occur in metal-mesh filters. 

The high-frequency edge of the observation band is 
defined by a 177 GHz (5.9 cm -1 ) multilayer LP filter 
mounted directly on the feed horn array at 280 mK. The 
low-frequency edge of the observation band is defined by 
1.33 mm diameter, 7.1 mm long cylindrical waveguides 
coupled to the feed horns described in Sec. |III C| Opti- 



cal transmittance within the cutoffs defined by the fil- 
ter and waveguide is further affected by the geometry 
of the bolometer absorption cavity. Measurements of the 
achieved band with a Fourier transform spectrometer are 
described in Sec. IWH1 
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280 mK 




Feed horn array 



TES bolometer array 



FIG. 4. Schematic diagram showing optical components in- 
side the APEX-SZ cryostat. Photons enter the cryostat via 
the Zotefoam window and pass through low-pass filters at 
60 K and 4 K, two HDPE lenses, and a band-defining filter at 
280 mK. Finally, smooth-walled conical horns couple photons 
to the bolometers. 

Each filter is clamped around its edge by an aluminum 
ring sandwiched with a beryllium-copper spring washer 
to compensate for differential thermal contraction of the 
aluminum and the polypropylene filter. The washer 
maintains thermal contact between the filter and ring 
at low temperatures. There is a 6-7 K temperature dif- 
ference between the front of the 4 K radiation shell and 
the center of the 255 GHz filter. Though the tempera- 
ture gradient is large, the filter does not add significant 
thermal load to the detectors, and the gradient itself is 
not a concern. 



C. Focal Plane Optics 

The APEX-SZ focal plane consists of 330 feed-horn 
coupled bolometers. Of these, 280 detectors are read out 
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Optical Element 


Diameter 
(mm) 


Shape 


Effective Focal 
Length (mm) 


Gaussian Beam 
Width (mm) 


Fractional 
Spillover 


Teriary mirror 


1440 


Off-axis paraboloid 


1970 


79 


0.001 


Folding flat mirror 


620 


Flat 


CO 


66 


0.001 


Reimaging mirror 


920 x 980 


Ellipsoid 


559 


58 


0.003 


4 K Lyot stop 


72 


Circular aperture 


CO 


45 


0.34 


Aperture lens 


120 


Double convex 


213 


47 


0.07 


Imaging lens 


240 


Double convex 


175 


36 


0.01 



TABLE II. APEX-SZ tertiary optics parameters. The optical elements are described in the text and in Fig. [2] The Gaussian 
beam width is the radius at which the power in the fundamental Gaussian mode for a single pixel's beam drops to e -2 of its 
peak power. The fractional spillover is fractional power vignetted by the optical element for a pixel on the edge of the focal 
plane array. 



(see Sec.JV]). The horn array mounted above the bolome- 
ter wafer is machined from a single aluminum block then 
gold plated. The horns are arranged in a hexagonal close- 
packed configuration with a 6.7 mm (1.4/A) spacing so 
that the entire focal plane is 133 mm in diameter. Fig- 
ure [5] shows a schematic diagram of the APEX-SZ inte- 
grating cavity; the design is similar to that of the Bolo- 
cam arrayF^ 

The detectors are thermally isolated, suspended struc- 
tures built upon a silicon wafer which is mounted on an 
invar plate. APEX-SZ has used two types of detectors 
which differ in the location of the reflective backshort. 
Type-1 detectors use the invar mount as a reflective back- 
short, and the wafer thickness, 450 /an, sets the back- 
short distance to be 3 A/4 in silicon. The type-2 detectors 
sandwich a metal film between two thin wafers to form 
a A/4 backshort 150 /um behind the bolometers. In both 
types, the space between the bolometers and backshort is 
mostly silicon. Section [VII| contains a detailed evaluation 
of the performance of the detectors. 

Each bolometer is centered behind a cylindrical waveg- 
uide in the horn array. Radiation not absorbed by the 
bolometer can leak radially to other detectors or be re- 
flected back through the horns. This leakage results in 
optical crosstalk between adjacent bolometers. We used 
a 3D-electromagnetic field simulation software package, 
High Frequency Structure Simulator (HFSS), 38 to sim- 
ulate and optimize the waveguide, cavity, and absorber 
geometries. Simulated results for the APEX-SZ config- 
urations are plotted in Fig. [6j The A/4 backshort has 
higher absorption and lower crosstalk across the obser- 
vation band. The band averaged crosstalk of the A/4 
geometry is 40% of that for the 3 A/4 geometry. 



IV. DETECTORS 

The APEX-SZ bolometers consist of an absorbing ele- 
ment coupled to a TES which is connected to a thermal 
reservoir by a weak thermal link. The TES is a thin su- 
perconducting film with a transition temperature higher 




FIG. 5. Schematic diagram of the APEX-SZ integrating cav- 
ity geometry for a single bolometer. The bolometers are 
lithographed on a silicon wafer and centered behind an ar- 
ray of cylindrical waveguides that couple to the smooth-wall 
horns. The horn/ waveguide array is positioned 400 /im from 
the wafer, and there is a 4 mm diameter, 250 (im deep relief at 
the end of the waveguide which reduces radiation loss along 
the wafer. The bolometer absorber is suspended on a silicon 
nitride membrane above a 20 /im vacuum gap etched in the 
silicon. The wafer is mounted on invar. In type-1 detectors, 
the invar mount serves as a 3A/4 backshort. In type-2 de- 
tectors, two thin wafers are joined with a metal film between 
them to serve as a A/4 backshort. The wafer and bolometer 
thicknesses are exaggerated for clarity. 



than the reservoir temperature. A constant bias volt- 
age is applied across the TES so the sum of electrical 
power dissipated in the TES plus optical power incident 
on the absorber heat the sensor to precisely the transi- 
tion temperature. Therefore, changes in the incident op- 
tical power produce an inverse change in electrical power. 
Since the voltage bias is held constant, the change in elec- 
trical power is measured as a change in current through 
the sensor. 
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FIG. 6. Results from HFSS simulation of the integrating cav- 
ity. For unit power input in the horn-coupled waveguide, 
the plot shows power absorbed by the spiderweb (solid lines) 
and power radiated radially in the wafer and air gap (dashed 
lines) as a function of frequency. The A/4 backshort geome- 
try (crosses) has better performance than the 3 A/4 backshort 
geometry (circles). 



The bolometer absorber is a 3 mm diameter gold spi- 
derweb, and the TES is located at its center (see Fig. [7| . 
The gold thickness of the absorber is set so the spider- 
web sheet resistance is 250 f2/sq. This resistance is a 
compromise between optimal absorption, which simula- 
tions suggest peaks at 500 £l/sq, an d fabrication con- 
straints, which require a minimum film thickness to main- 
tain metal continuity across the absorber. The central 
location of the TES and the low heat capacity of the 
micromesh spiderweb absorber result in a fast optical re- 
sponse. The properties of the bolometers are summarized 
in Table Cni 



Parameter 


Type-1 


Type-2 


Rn (n) 


1.9 


1.1 


T c (mK) 


465 


570 


G (pW/K) 


153 


86 


T opt (ms) 


11 


20 


Responsivity (xl0 5 A/W) 


1.7 


3.0 



TABLE III. Typical TES Bolometer Parameters, 
are defined in the text of Sec. II VI 



The terms 



The absorber, TES, and leads to bias the TES are de- 
posited onto a 1 /im thick low-stress silicon nitride (LSN) 
spiderweb membrane. The LSN spiderweb is suspended 
by eight LSN legs above a ~ 20 ^m vacuum gap etched 
in the silicon. The gap is made by a dry XeF2 gas etch, 
and thermally isolates the absorber from the silicon sub- 
strate. The thermal conductance of the bolometer is set 



by a gold thermal link that runs from the TES to the 
thermal reservoir, parallel to the bias leads. The sup- 
port legs are an additional, low conductivity link to the 
thermal reservoir. A fraction of the incident power ab- 
sorbed on the spiderweb reaches the reservoir via the legs 
and does not heat the sensor as much as power reaching 
the reservoir via the gold link. A finite element analysis 
of the thermal structure of the bolometer and the ex- 
pected illumination pattern gives a calculated bolometer 
absorption efficiency, r\boio, of 0.84 wh ich low ers the effec- 
tive efficiency of the bolometer (Sec. VII F). The thick- 
ness of the gold link is tuned so that the total thermal 
conductance to the reservoir, including the LSN legs, is 
G = 153 pW/K (for type-1, see Table |lu). Here, G is 
defined as 



G = 



P 



bias 



P 



opt 



T r -T< 







(1) 



where T c is the superconducting transition temperature 
is, To is the reservoir temperature (280 mK in APEX- 
SZ), and Puas and P op t are, respectively, the total elec- 
trical and optical power dissipated in the bolometer. 

For the detector, T c and G must be carefully chosen to 
avoid saturation while maximizing sensitivity. Following 
equation (JlJ, T c and G must be large enough to dissi- 
pate the expected optical power with enough margin to 
allow for loading variations and sufficient bias power. We 
typically use 



G(T C - T ) « 2P t 



opt ■ 



(2) 



The minimum thermal carrier noise for a bolometer 
where electrons are the dominant thermal carrier is given 
by the conditional 



T B = 2.1 T . 



(3) 



The minimum thermal carrier noise gives the maximum 
contribution by photon noise to the total noise of the 
bolometer. With To fixed by our choice of cooler, equa- 
tions ([2]) and ([3]) determine the appropriate T c and G 
for an estimated P op t- In practice, achieving a specific 
transition temperature is difficult and optical power is 
not precisely known at the detector fabrication stage. As 
discussed in Sec. 



VII H 



the type-2 detector parameters 
result in increased sensitivity relative to the type-1 de- 
tectors. 

The TES is an Al-Ti bilayer with layer thicknesses 
tuned for a T c of 465 mK and geometry chosen for a 
normal resistance, R n , of 1.9 f2(for type-1). The TES is 
voltage biased and operated with strong electrothermal 
feedback (ETF). The ETF maintains the TES in the su- 
perconducting transition and increases the range of the 
linear response.^ The loop gain of the ETF is given by 



£ = 



Pbias®- 

GTr 1 



where Pu as is the bias power and a = dlog(R)/dlog(T) 
is a measure of the sharpness of the superconducting 



FIG. 7. (a) Three spiderweb bolometers on a 55-element sub-array, (b) A close-up of the center of the spiderweb shows the 
elements of the bolometer. The dark areas along the absorber legs and inside the gold ring are regions in which the silicon has 
been etched away from the spiderweb. 



transition. Note that G refers to the instantaneous ther- 
mal conductivity of the bolometer, 5P/5T, which de- 
pends on the material properties of the link. For a detec- 
tor biased with a constant voltage, the current respon- 
sivity is determined by the voltage bias 



SP V Mas C + 1 V bias ' 

For the sinusoidal bias system used in APEX-SZ, the 
responsivity is modified by a factor y/2, but the NEP for 
constant and alternating-biased systems are identical in 
the ideal case as discussed in Sec. |Vl 

The ETF also speeds up the intrinsic thermal response 
time of the bolometer, tq = C/G where C is the heat 
capacity of the detector, so the thermal time constant is 

TO 

Tth -£Ti- 

This thermal time constant sets the TES response time. 
In the APEX-SZ bolometers, the optical response time 
is limited by the thermalization time of the spiderweb 
absorber, T opt > r t h- The T opt of the type-1 and type-2 
detectors is 11 ms and 20 ms, respectively. The long time 
constant of the type-2 bolometers is likely due to a pro- 
cessing anomaly, high-heat capacity, optically-inactive 
residue on the spiderweb. We must take the long time 
constant into account in data analysis, but the sensitivity 
is not affected. 

For stable voltage bias, we require that the TES ther- 
mal time constant be larger than the time constant of 
the bias circuit. The stability criteria given by Ref. l4"Tlis 

Tth/5.8 > Tbias- 



The bias circuit time constant is dominated by the LCR 
tank circuit of the readout described in Sec. [V] and is a 
function of the TES resistance, Tbi as = L/R. The intrin- 
sic TES thermal time constant is < 100 [is, too fast for 
stable operation at high ETF loop gain and much faster 
than the optical time constant from thermalization of 
the spiderweb. Therefore, we slow the thermal response 
of the bolometer by adding additional heat capacity in 
the form of a 200 /im wide, 3 /im thick gold ring coupled 
to the TES. The ring geometry is chosen to allow stable 
bolometer operation at loop gains of C = 10-100 while 
still keeping the bolometer thermal time constant much 
less than the bolometer optical time constant. 



The APEX-SZ detector array is assembled from six 
identical 55-element, triangular sub-arrays (Fig. [8]). Each 
sub-array is fabricated on a single 4 inch silicon wafer 
with standard optical photolithographic techniques at 
the Berkeley Microfabrication Laboratory. The triangu- 
lar sub-arrays are mounted on matching invar triangles 
with BeCu spring clips at the corners. Invar is used to 
match the thermal contraction of the wafer. Thermal 
contact between the wafer and invar is made by a thin 
film of Apiezon-N greaseP^I The invar is itself attached 
to an aluminum mount which supports a multiplexer cir- 
cuit board, and the bolometer sub-array is wirebonded to 
the board. In the receiver, the triangular sub-assemblies 
form a single planar hexagonal array with 330 bolome- 
ters. The modular design allows construction of large 
focal plane arrays and easy replacement of sub-arrays. 
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25 mm 




FIG. 8. (a) A fabricated bolometer sub-array in its holder. Each 55-element sub-array is mounted on an invar triangle and 
wirebonded to a multiplexer circuit board. Multiplexer circuit elements (tuned capacitors and inductors, described in Sec. [V| 
are housed under the aluminum guard. A single micro-D connector provides the connection to the SQUIDs. (b) The TES 
bolometer array. Six identical sub-arrays are assembled into single planar array which has 330 bolometers and is 133 mm in 
diameter. Each bolometer has two leads which run from the spiderweb to bonding pads at the edge of each triangle. These are 
visible as light traces between the bolometers. 



V. FREQUENCY MULTIPLEXED READOUT 
ELECTRONICS 

The APEX-SZ TES bolometers are read out with a 
SQUID amplifier frequency-domain multiplexing (FDM) 
system 43 which allows a number of detectors to be read 
out by a single 4 K SQUID amplifier connected through 
a single pair of wires. Readout multiplexing greatly re- 
duces the thermal load on the low-temperature stage, 
the complexity of cold wiring, and the system cost. This 
technique enables the use of large detector arrays. The 
APEX-SZ system uses seven detectors read out by one 
SQUID amplifier in each multiplexer module. 

The basic components of a readout system module are 
shown in Fig. [9] The bolometers are biased with alternat- 
ing voltages at carrier frequencies (0.3-1 MHz) that are 
much higher than the bolometer thermal bandwidth, so 
the bias deposits a power on the sensor that depends only 
on its resistance. Each bolometer in the module is biased 
at a different frequency. Changes in absorbed radiation 
change the resistance of a given bolometer and modulate 
the current. This amplitude modulation translates the 
absorbed signal spectrum to sidebands centered around 
the carrier bias frequency. Since the currents from the 
different bolometers are separated in frequency, they can 
be combined and transmitted to the SQUID amplifier 
with a single wire. Furthermore, the bolometers are con- 
nected through series LC circuits, each of which is tuned 
to the appropriate bias frequency, so all bias frequen- 
cies can also be fed through a single line. As a result, 
each multiplexer module requires only a single pair of 
wires from the low-temperature stage to the 4 K stage 
to read out all the bolometers. The comb of amplitude 
modulated carriers is fed to a bank of demodulators that 
mix the signals down to base-band. The signals are then 



filtered and digitized. All outputs in the array are sam- 
pled synchronously at 1 kHz, a much higher frequency 
than the signal bandwidth. The data streams are passed 
though a digital low-pass filter, downsampled to 100 Hz, 
and compressed before being written to disk by the read- 
out computer. 

High impedance NTD bolometers are susceptible to 
spurious signals induced by microphonic vibrations^ of 
the bolometer wiring. However, TES bolometers are low 
impedance devices which are relatively insensitive to vi- 
brations. Using an alternating voltage bias further re- 
duces the sensitivity to microphonics. Since the detec- 
tors only respond to voltages near their bias frequency, 
the FDM system is insensitive to microphonic excitation 
of the bolometer wiring which occurs at much lower fre- 
quencies. Thus, no special effort to restrain the bolome- 
ter wiring is required for this system. 

The bandwidth of the LC circuits is chosen to accom- 
modate the signal sidebands and attenuate the crosstalk 
between neighboring channels to an acceptable level (< 
1%). The bandwidth is the same for all channels and 
is determined by the inductance and the bolometer re- 
sistance, Rboiol '(2-7T-L). The high-Q LC filters are com- 
posed of photolithographically constructed 16 /iH spiral 
inductors and commercial negative-positive-zero (NPO) 
ceramic chip capacitors that are mounted on a printed 
circuit board adjacent to the bolometer sub-array. The 
resulting filter is compact and comparable in area to a 
single bolometer. In addition to defining the carrier fre- 
quency for each channel, the tuned circuits also limit the 
bandwidth of the bolometer Johnson noise, which would 
otherwise contribute to the noise in all other channels of 
the module. 

The SQUID amplifier system has been carefully de- 
signed to achieve the required low input impedance, low 
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FIG. 9. Schematic diagram showing a single module of the frequency multiplexed SQUID readout system. Sinusoidal bias 
voltages at different frequencies are summed and applied to a group of bolometers via a single line. LC filters in series with 
the bolometers isolate an individual bias for each bolometer. The modulated signals from the bolometers are summed at the 
input of a SQUID array. 



noise, and high dynamic range across the MHz car- 
rier bandwidth. Low input impedance is necessary be- 
cause voltage biasing requires that the impedance of the 
SQUID preampliflier be much smaller than the bolometer 
resistance. Shunt feedback is applied to further reduce 
the input impedance of the preamplifier. To achieve the 
necessary transimpedance and dynamic range, an inte- 
grated array of 100 SQUIDs with inputs and outputs con- 
nected in series^ is employed. The SQUID arrays were 
designed and fabricated at NISTP The array allows a 
smaller total input inductance and increased signal-to- 
noise relative to a single SQUID. A feedback loop gain 
of about 20 is achieved by feeding the SQUID output to 
a room-temperature amplifier with low noise and a suf- 
ficiently large bandwidth to maintain stability over the 
full range of bias frequencies. 

The bolometer signal information appears in the side- 
bands of the carrier frequency while the amplitude of 
the carrier contains no useful information. Thus, the dy- 
namic range requirements of the SQUID system can be 
reduced by nulling the carrier without losing any signal 
information. At each carrier frequency, we add a constant 
amplitude nulling current which is adjusted to be out of 
phase by tt at the SQUID input. The resulting nulling 
of the carrier reduces the total signal amplitude through 
the SQUID amplifier system by more than a factor of 10. 
The bias and nulling signals take different paths through 
the cryogenic system and so they incur different phase 
shifts due to inductive and capacitive strays such as the 
wiring inductance and the reactance of the LC filters on 



neighboring channels. 

To tune the phase of the miller waveforms efficiently 
it is necessary to measure the magnitude of the resid- 
ual carrier. This measurement is achieved by using an 
eighth multiplexer channel from each module to measure 
the magnitude of the residual bias for the other channels 
during the nulling process. This strategy improves the 
speed and accuracy of the carrier nulling at the expense 
of a channel which would otherwise be used to readout an 
additional detector. Each multiplexer module has eight 
channels, seven used for bolometer readout and one used 
for nulling. 

The seven bolometers per module are read out by a 
single SQUID amplifier. For each triangular sub-array 
of bolometers, seven groups of seven bolometers are each 
fed to a single SQUID card which holds the seven SQUID 
arrays. Thus, the readout requires fourteen wires for each 
sub-array. Note that there happen to be seven bolome- 
ters for each SQUID array and seven readout SQUID 
arrays on each card, but these two factors are indepen- 
dent. The wiring must have low thermal conductance 
to minimize the load on the millikelvin stages, and low 
impedance to minimize shifting of the LC resonance fre- 
quency. For this purpose, we use a hybrid wiring shown 
in Fig. |10| low inductance copper traces coated with tin- 
lead solder on flexible kapton in series with low thermal 
conductance NbTi twisted-pairs. The traces are 1.8 mm 
wide, deposited on both sides of a 18 mm wide kapton 
ribbon, thereby limiting the inductance to 0.05 nH/mm. 
The NbTi wires have a CuNi coating for soldering, are 
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woven into a cabled and have 1 nH/mm inductance. 
Both parts of the hybrid cable are superconducting. Us- 
ing 510 mm of flexible-ribbon cable with 100 mm of 
twisted NbTi cable provides an acceptable compromise 
between inductance, ~ 150 nH total, and thermal load 



(Sec. |VID| ). 

In the APEX-SZ configuration, we use six SQUID 
cards each with seven SQUID arrays to read out bolome- 
ters. Each SQUID card also has an eighth SQUID array 
which is not wired to bolometers and can be used for 
readout diagnostics. The SQUIDs are extremely sensitive 
to magnetic fields and must be carefully shielded. Each 
SQUID is mounted on a 9 mm square pad of supercon- 



ducting niobium foil (Fig. 10 ), which serves to trap resid 



ual flux during cool down and attenuate time-varying 
magnetic fields. Further, the SQUID card is housed in- 
side its own CryoperrrpS magnetic sheath. Cryoperm is 
an alloy which has high permeability at cryogenic tem- 
peratures. The relative permeability at 4 K is 70, 000 
for 0.4 A/m static fields, comparable to the permeabil- 
ity of mu-metal at room temperature. The Cryoperm 
shield attenuates both time-varying and spatially vary- 
ing fields. We measure a reduction of field variations 
by a factor 2.5xl0 4 with this two-fold shielding scheme. 
At this shielding level, SQUID movement through the 
Earth's magnetic field is detectable at signal-to-noise ra- 
tio of near unity. However, such low spatial frequencies 
are filtered by the analysis pipeline. 

As noted above, the high gain-bandwidth product of 
the feedback loop requires short connections between the 
4 K SQUID arrays and room temperature op-amps to 
minimize phase shifts and maintain stability of the neg- 
ative feedback loop at megahertz frequencies and high 
loop gains. The amplifiers arc mounted on "SQUID 
controller" circuit boards which are attached directly 
to the side of the receiver cryostat wall in an RF-tight 
aluminum box. There is a one-to-one correspondence 
between SQUID boards and SQUID controllers; each 
SQUID controller board houses the op-amps and con- 
trol logic to provide bias currents for the eight SQUID- 
arrays on each SQUID board. To couple the 4 K SQUID 
boards to the room temperature SQUID controllers, we 
use a custom manufactured^ wire harness consisting of 
518 120 mm length manganin wires in woven Nomex. Up 
to seven SQUID cards can be plugged into a single wiring 
harness. 

The comb of sky-signal modulated carriers is passed 
from the SQUID controller to a bank of demodulators on 
one of twenty 16-channel oscillator-demodulator boards. 
Direct Digital Synthesizers (DDSs) provide carrier and 
nulling signals with precise frequency and amplitude con- 
trol and very low sideband noise. The same DDS that 
generates the bolometer bias also provides a phase and 
frequency locked square-wave reference for the corre- 
sponding demodulator. In the absence of phase shifts in 
the cryostat, the demodulator output would correspond 
to the in-phase /-component of the carrier signal. How- 
ever, the wiring strays result in a small offset in phase 



between the demodulator and carrier that increases with 
frequency. The lack of phase adjustment in the demod- 
ulator means this phase contributes to a slight degrada- 
tion of noise performance at high carrier frequency (see 
Sec. VII H |. The demodulator circuits use a sampling de- 



modulator followed by an 8-pole, low-pass anti-aliasing 
filter, and a 14-bit analog-to-digital converter (ADC). 
A field programmable gate array (FPGA) assembles the 
data from all channels on a given board and streams it to 
the data acquisition computer. Each demodulator board 
can process the outputs from two SQUID amplifiers, i.e. 
two groups of eight multiplexer channels. The readout 
crate can accommodate twenty demodulator boards with 
a total of 320 channels. However, since only seven chan- 
nels in each group are used to read bolometers, the max- 
imum number of bolometer channels is 280. There are 
seven channels for each group of eight bolometers. 



VI. CRYOGENICS 
A. Cryostat 

The APEX-SZ cryostat houses cold filters and lenses, 
the detector array, and SQUID electronics. The focal 
plane is cooled by a 3-stage helium sorption refrigerator 

sum- 



backed by a mechanical pulse-tube cooler; Tabic IV : 
marizes the available cooling power at various tempera- 
ture stages. The mechanical cooler eliminates the need 
for open reservoirs of liquid cryogens, greatly simplifying 
the design and construction of the cryostat while reduc- 
ing the cost and difficulty of remote observations. The 
receiver is essentially downward-looking, further simpli- 
fying the cryostat design. The optical and pulse-tube 
axes are parallel while maintaining the pulse-tube cooler 
within 30° of vertical orientation as required for its op- 
eration. 



Stage 


Cooling Power (W) 


60 K 


40 


4 K 


1 


350 mK 


6.0xl0 -5 


270 mK 


1.5xl0~ 6 



TABLE IV. Cooling power provided by the pulse-tube cooler 
and the 3-stage helium sorption refrigerator at different tem- 
perature stages of the APEX-SZ cryostat. 



Cutaway drawings of the inverted cryostat are shown 
in Fig. [TTJ The cryostat includes 4 K and 60 K radia- 
tion shields supported by G-10 fiberglass struts inside the 
300 K vacuum shell. The 60 K radiation shield is gold 
plated to reduce its emissivity and decrease the radiation 
load on both the 4 K and 60 K shields. Where feasible, we 
wrap cylindrical portions of both shells with aluminized 
mylar to further reduce the radiative load. The optical 
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FIG. 10. (a) Two SQUID cards with 4 K to 280 mK wiring and a Cryoperm sheath. Traces on flexible kapton ribbon, soldered 
directly to the SQUID card, are in series with woven NbTi cable. In APEX-SZ we use 510 and 360 mm lengths of flexible-ribbon 
cable. The NbTi cable has a copper block epoxied at its midpoint for heatsinking at 370 mK. The male micro-D connector at 
the end of the cable mates with the female connector on the multiplexer board shown in Fig. [8] (b) A detailed image of the 
SQUID card with magnetic shielding. The SQUID card has 8 SQUID arrays, each of which is mounted on a square pad of 
niobium foil. The card slides into the Cryoperm sheath shown above the card, leaving the gold contacts exposed for inserting 
into a Peripheral Component Interconnect (PCI) connector. The left side of the card has 8 bias chip resistors. The upper edge 
of the photograph shows part of a SQUID card before the SQUID arrays have been attached, (c) A close-up image of the 
100-element series array SQUID chip. 



section can be removed for access to the focal plane and 
SQUID electronics. The section of the receiver which 
houses the wiring modules was designed to be modular, 
allowing aspects of the cryostat design to be used for 
future experiments. 

Because the APEX telescope hosts a number of in- 
struments, the uncertain RF environment of the cabin 
was a major concern. The lower section of the cryostat 
(as oriented in Fig. 11) is an RF sealed enclosure. All 



vacuum o-ring seals also include an additional RF seal 
of elastofoamPS SQUID controller cards are housed in 
an RF tight box mounted on the cryostat, with boards 
plugging directly into the hermetic feedthroughs on the 
cryostat shell. Commercial pi-filter connectors on the 
outer face of the box provide feedthroughs for the read- 
out wiring. Inside the cryostat, aluminized mylar sheets 
create electrical shorts between the 300 K, 60 K and 4 K 
shields. An additional sheet between 4 K and 280 mK 
completes the shielded compartment. The only perfora- 
tions are the feed horn array, hollow screws to connect 
vacuum spaces, and the electrical feedthroughs on the 
room-temperature electronics box. The feed horns are 
coupled to waveguides with a 130 GHz cutoff frequency, 



the holes in the hollow screws have a high length to diam- 
eter ratio, and the electrical feedthroughs arc protected 
by pi-filters. The thickness of the aluminum layer on the 
mylar is a compromise between electrical and thermal re- 
sistance, with better shielding also increasing the thermal 
load on low temperature stages. We take advantage of 
the larger cooling power available at higher temperatures 
and tailor the aluminum thickness for each temperature 
stage: 400 nm for the 300-60 K temperature gap, 100 nm 
for 60-4 K, and 50 nm for 4-0.280 K. 



Heat straps from the cryostat cooling elements — the 
pulse-tube cooler heads and sorption refrigerator mil- 
likelvin stages — to the cryostat are all made from lengths 
of 12 AWGP braided oxygen-free high conductivity 
(OFHC) copper wire melted into solid copper blocks at 
each end by tungsten inert gas (TIG) welding. Braided 
wire was chosen to decouple the vibrations of the pulse- 
tube from the receiver and allow differential thermal con- 
traction between the cryostat components. The solid 
copper ends give high thermal conductance interfaces. 
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FIG. 11. Drawings showing two orthogonal sections of the APEX-SZ cryostat. Many components are labeled. The cryostat is 
shown looking upward in the orientation used in the laboratory for work on the focal plane. The optical section (upper half) 
can be removed for access to the focal plane and readout components. The aluminized mylar sheet (red) acts an RF shield and 
separates the optical and readout sections of the cryostat. Optical components are shown in light blue, 280 mK components in 
orange, 380 mK in yellow, and the 4 K SQUID card housing and readout wiring module in green. When mounted in APEX, 
the cryostat is downward-looking, with the window pointed at a mirror on the floor of the cabin. 



B. Pulse- Tube Cooler 

APEX-SZ uses a commercial pulse-tube refrigerator, 
model PT410-OP3 manufactured by Cryomech, Inc.f^to 
provide cooling at 4 K and 60 K. In a pulse-tube cooler, 
high pressure He gas from a compressor is passed through 
a high heat capacity regenerator and then allowed to ex- 
pand while in contact with a cold head, absorbing heat 
from the head. A motor outside the cryostat drives a 
valve which controls the pressure cycle of the pulse-tube. 
There are no moving parts or valves in the cold section. 
This feature makes the pulse-tube cooler lower vibration, 
more robust, and easier to maintain than other mechan- 
ical coolers. 

The pulse-tube cooler has a number of advantages over 
liquid cryogens for remote observations, but it can add 
noise to measurements via several mechanisms: electri- 
cal noise, mechanical vibration, and temperature oscilla- 
tions. The compressor and valve motor are electrically 
noisy and can cause oscillating voltages on the cryostat 
shells. The cooler used in APEX-SZ has its valve mo- 
tor separated from the pulse-tube by a 500 mm length of 
flexible tubing including a short section of HDPE which 
electrically isolates the compressor and valve motor from 
the cryostat. In addition to electrically isolating the 
valve motor, we drive it with a low-noise linear stepper 
motoric Valve motor movement and gas flow propagate 



vibrations along the pulse-tube and cause the cold head 
to vibrate. The pulse-tube is mechanically isolated by 
mounting it on an elastomer vibration damper. 54 Flexi- 
ble copper braids isolate components inside the cryostat 
from vibrations on the cold heads. The temperature of 
the cold head oscillates with the pressure cycle. The bare 
pulse-tube cooler head has 200 mK temperature oscilla- 
tions at the 1.4 Hz cycle frequency but these are atten- 
uated by the large heat capacity of the cryostat. The 
resulting oscillations at different cryostat stages are de- 
tailed in Sec. IVTEl 



An additional concern when using the pulse-tube 
cooler rather than liquid cryogens is that the cooling per- 
formance varies with the orientation of the pulse-tube 
relative to gravity. The best performance is achieved 
with the pulse-tube pointing downward. The PT410 loses 
about 10% of its cooling power when tilted 30° from verti- 
cal, with more severe losses at larger angles. The tertiary 
optics, cryostat and mount were designed to hold the 
cryostat and pulse-tube at 30° when pointed at zenith. 
Thus, for observations between zenith and 30° elevation, 
the pulse-tube is tilted at less than 30° relative to verti- 
cal. 
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C. Helium Sorption Refrigerator 

We use a custom built 3-stage helium sorption refrig- 
erator from Chase Research^ to cool the bolometers to 
280 mK. It is sometimes referred to as a "He- 10" refrig- 
erator because it includes one 4 He stage and two 3 He 
stages. The design includes sintered copper in the he- 
lium reservoirs to provide good thermal contact between 
the fluid helium and the refrigerator cold heads over a 
range of orientations. Many previous bolometric CMB 
experiments such as the Arcminute Cosmology Bolome- 
ter Array Receiver ( ACB AR; Ref. [55] and Bolocam have 
used similar refrigerators backed by open reservoirs of 
liquid cryogens. 

A nominal Hc-10 refrigerator cycle cools the APEX-SZ 
focal plane to 280 mK for more than 24 h. The 4 He is cy- 
cled only to condense 3 He into the "Interhead" and "Ul- 
trahead" reservoirs which cool to 350 mK and 270 mK, 
respectively. At these temperatures (and ignoring para- 
sitic loads), the Interhead has 60 fiW of cooling power, 
absorbing 5 J during the cycle. It serves to buffer the Ul- 
trahead, which has 1.5 fiW of cooling power, absorbing 
0.1 J during the cycle. Heat from the wiring, RF shield- 
ing and focal plane support structure is intercepted at 
the Interhead. 

The efficiency of the refrigerator cycle as a whole de- 
pends critically on the temperature of the condensation 
point of the 4 He stage. It is therefore thermally con- 
nected directly to the 4 K head of the pulse-tube via 
braided cooper wire as described in Sec. |VI A| 

The He- 10 refrigerator is cycled by applying currents 
to six resistors. Three operate the charcoal pumps, and 
three operate gas-gap heat switches between the pumps 
and 4 K head. The switches are closed by using a resistor 
to heat a small charcoal pump and release helium into the 
space between the hot and cold surfaces of the switch. An 
electronics box with computer control logs temperatures 
(Sec. VI E| and allows remote, automatic cycling of the 
refrigerator in 2.5 h. At the end of the cycle, we bias the 
bolometers for observation while they are in the normal 
state and then allow the focal plane to cool below T c with 
the bolometers biased. 



D. Millikelvin Stage 

The support structure for the focal plane provides a 
rigid mechanical support while maintaining the thermal 
isolation required to achieve millikelvin temperatures. As 
shown in Fig. [IT] the focal plane support consists of three 
stages supported from the 4 K stage. The structure pro- 
vides temperature intercepts with each stage thermally 
anchored to a different point of the He-10 refrigerator: 
the ~ 3 K stage is thermally connected to the "Heat Ex- 
changer" of the refrigerator (a point which equilibrates 
between the 4 K stage and Interhead temperatures); the 
~ 370 mK stage is thermally connected to the Interhead, 
and the ~ 280 mK stage is thermally connected to the 



refrigerator Ultrahead. The focal plane and LC filter 
boards are bolted to the ~ 280 mK stage. Wires to the 
focal plane and the aluminized mylar RF shield are ther- 
mally anchored at the ~ 370 mK stage. The NbTi section 
of low-inductance hybrid wiring (described in Sec. |v| is 
thermally anchored at the ~ 370 mK stage resulting in 
loads of 16 fiW on the Interhead and 0.03 fiW on the 
Ultrahead. 

The stages are made from gold-plated aluminum, and 
each level is supported from the one below it by six 
VespeP^ legs. We use Vespel SP-1 for supporting the 
~ 3 K and ~ 370 mK stages. The load from the sup- 
ports on the Interhead is 19 /zW. We use graphite loaded 
Vespel SP-22 to support the ~ 280 mK stage due to its 
lower thermal conductance compared to unloaded Vespel 
at millikelvin temperatures.^ The load on the Ultrahead 
from the supports is 0.3 /iW. The full thermal loading on 
the two heads is summarized in Tabic IVl 

Interhead (jiW) Ultrahead (/iW) 



Support structure 
RF Shielding 
Wiring 
Radiation 
Bolometers 

Total 



19 
15 
16 
0.004 



50 



0.32 
0.15 
0.03 
0.19 
0.01 

0.70 



TABLE V. Power Dissipation on mK Stages 



E. Thermometry 

APEX-SZ uses commercial thermometers^ read out 
by the data acquisition computer to monitor temper- 
atures inside the cryostat. The refrigerator millikelvin 
heads and their corresponding stages are monitored with 
Cernox RTD thermometers, accurate from 300 K to 
below 280 mK and read with a AC bridge. Silicon 
diode thermometers monitor refrigerator pumps and heat 
switches as well as both the 60 K and 4 K pulse-tube 
heads and corresponding stages. The temperatures re- 
turned from the thermometers as well as the states of 
the six heaters on the sorption refrigerator are all contin- 
uously monitored via an electronics box in the Cassegrain 
cabin similar to that used in ACBAR. 

We do not use active temperature stabilization at 4 K 
or on the focal plane. Slow temperature drifts of the 
focal plane or the optics have an observable effect on 
the bolometer signal, but are well below the signal band. 
The pulse-tube cooler has 1.4 Hz temperature oscilla- 
tions which cause 2 mK fluctuations at the 4 K cryostat 
stage. The temperature fluctuations at the Ultrahead 
are < 20 nK. These fluctuations induce a ~ 2 aW sig- 
nal in the bolometer which is well below the detector 
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noise level. Temperature fluctuations in the cold optics 
cause variable loading on the bolometers, but these are 
damped by the thermal resistance along the radiation 
shells and the heat capacity of the optics. The Lyot stop 
has temperature fluctuations < 3 /iK, roughly equivalent 
to a 3 /iK CMB temperature fluctuation which is well 
below the detector noise. Because the scan strategy is 
asynchronous with pulse-tube fluctuations, any resulting 
signal would appear as an insignificant additional noise 
in the final maps. 



channels. After data processing, there are typically 120- 
140 type-1 channels which contribute to a cluster map. 

On the type-2 sub-array installed in 2009 and used 
in this analysis, 70% of the bolometers can be stably 
biased in the transition. After data processing, there are 
typically twenty type-2 detectors which contribute to a 
map. 

B. Beams 



VII. SYSTEM PERFORMANCE 

The APEX-SZ receiver was first installed at the APEX 
telescope in December 2005 for an engineering run with a 
single 55-element sub-arrayP^ We redeployed the system 
with a full 280 channel array of type-1 detectors in early 
2007. In 2009, we replaced one type-1 sub-array with a 
type-2 sub-array. 

As a millimeter-wave instrument that is nominally less 
affected by atmospheric water vapor than shorter wave- 
length receivers, APEX-SZ has observed in a range of 
weather conditions. In bad weather, detector sensitiv- 
ity can degrade by up to a factor two, and aggressive 
timestream data filtering is required to eliminate large 
spatial-scale noise due to temporally varying water vapor 
fluctuations. A full characterization of the instrument 
and the Atacama site in various weather conditions is 
beyond the scope of this paper. Therefore, as an accu- 
rate assessment of the receiver's capabilities, we analyze 
data taken during excellent weather at the height of the 
Austral winter. In this section, we use data from an 
August 2007 observing period to characterize the opti- 
cal performance of the system and the sensitivity of the 
type-1 detectors. The type-2 detectors are characterized 
using data from April 2009. 



A. Detector Yield 

Typically, 170-180 of the 280 type-1 channels in the re- 
ceiver were optically responsive. One of the sub-arrays, 
fabricated earlier than the others, was unstable and had 
no live bolometers (as indicated by the wedge-shaped 
void in Fig. 12 1. Two SQUIDs had anomalously high 
noise, which resulted in the loss of 14 channels. In the 



other five sub-arrays, fabricated with a more mature pro- 
cess, 90% of the TESs could be biased in their transition 
with strong ETF. 

The gold layer on the spider-web absorber is as thin as 
the manufacturing process allows to match the charac- 
teristic wave impedance in the cavity. It is likely several 
bolometers have incomplete gold coverage which reduces 
optical absorption. About 85% of the bolometers which 
can be stably biased are optically sensitive enough to 
map bright calibration sources. However, noise filtering 
results in removing additional noisy and low sensitivity 



The APEX-SZ beams are measured with a raster scan 
of a planet. The rasters are 0.5° x 0.5° with azimuthal 
sweeps separated by 18" in elevation. This scan pattern 
is wide enough to sample the full array and dense enough 
to ensure that each beam is fully sampled. 



1. Gaussian Fit 

For each channel, timestream data are binned into an 
individual beam map. The planet source is frequently 
extended in our beam, therefore, each map is fit to the 
convolution product of a 2-D Gaussian and a planet-sized 
disk. The resulting best-fit 2-D Gaussian determines the 
size, ellipticity and pointing offset for each beam. 

Figures [12] and [13] show the result from a typical Mars 
scan. During this scan, the angular size of Mars was 8". 
The median FWHM beamwidth is 58 ± 6" with a range 
of 45-75" and a median axial ratio of 1.17. The median 
beam area is 1.0 ± 0.1 amin 2 . The measured FWHM 
beamwidths and axial ratios agree quite well with those 
expected from the optics design. A sample of nine beams 
across the array modeled using ZEMAX physical optics 
propagation have median FWHM beamwidth of 57" with 
range of 47-75" and median axial ratio of 1.09. The beam 
size, ellipticity and orientation change slightly across the 
array, with ellipticity increasing toward the edges. 

The size and orientation of the outer beams is sensi- 
tive to the focus position of the secondary mirror. The 
best focus position is that which minimizes the median 
beam FWHM across the array. This position is slightly 
offset from the best focus position for a channel near the 
array center. We typically measure this offset distance 
at the beginning of an observation period by making a 
set of beam maps at various focus positions. For subse- 
quent days, we find the best focus position for a single 
central pixel and apply the secondary offset. The offset 
is remarkably consistent between observing periods. 



2. Beam Profile 

In addition to a central Gaussian lobe, each beam 
has an extended non-Gaussian profile. Mapping the 
extended profile for each beam would require a pro- 
hibitively long scan, over 50 h of integration on a calibra- 
tion source. Therefore, we create a coadded map of the 



16 



CD 

tn 
o 



LU 



600 
400 
200 


-200 
-400 
-600 



®A l© n 



% 

©©(+ 



®©©© @ ® © 

«£® @ ©\ 

%%®% ©% @ © 

O 



-600 -400 -200 200 400 

A AZ (arcsec) 



600 800 



FIG. 12. Elliptical Gaussian beam fits of 177 optically live 
channels. Ellipses represent best fit 2-D Gaussian FWHM of 
individual channel maps; pluses mark the center of each beam. 
The circle in the lower left corner is 60". The ellipticity and 
orientation of the beams change slightly across the array. The 
median beamwidth is 58±6" with a median axial ratio of 1.17. 



profile. 

The measured mean beam profiles for both detector 
types have a central lobe that is best fit by a 59" Gaus- 
sian. The profiles differ in the level of the near sidelobes, 
— 14 and — 15 dB for the type-1 and type-2 detectors, re- 
spectively. The sidelobes increase the beam solid angle 
compared to the best fit Gaussian by 30%. The effective 
beam area relative to the Gaussian is factored into the 



calibration of the instrument (Sec. VII D ) 



We compare the beam profiles of both detector types 



to the theoretical prediction in Fig. 14 Theoretical beam 
shapes for different positions in the array are calculated 
using zemax. The theoretical profile in Fig. [14] is calcu- 
lated from a weighted composite of the predicted beams 
from several positions in the array. The predicted beam 
has sidelobes due to truncation at the Lyot stop. The 
measured beams show increased near-sidelobe power over 
that predicted which is consistent with optical crosstalk 
between adjacent bolometers. The sidelobe level of the 
type-1 (type-2) profile indicates 2.5% (< 1%) crosstalk 
between adjacent bolometers. This level of crosstalk and 
the factor 2.5 difference in crosstalk is consistent with 
cavity simulations (Fig. [6]). 
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FIG. 13. The lines show the mean Gaussian beam area (black) 
and axial ratio (blue) versus radial offset from the focal plane 
center. Error bars show the standard deviation of each value 
for the pixels in each radial bin. The beam area and axial 
ratio increase by 17% from the center to edge of the array, 
but the variation in beam area increases by a factor five. 



average beam by combining individual beam maps with 
minimum variance weighting, using the sample variance 
of the timestream data. Though the individual beams 
are elliptical and their orientation varies, the composite 
map is very nearly radially symmetric. We take the ra- 
dial average to derive a single high signal-to-noise beam 
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FIG. 14. Comparison of measured versus predicted average 
beam pattern of type-1 and type-2 detectors. The individual 
channel maps from a planet observation are coadded and radi- 
ally averaged to derive an average beam profile for the type-1 
(black solid line) and type-2 (red solid line) detectors. The 
average predicted beam, also radially averaged (dashed lines), 
shows near sidelobes due to truncation at the Lyot stop. The 
measured beam central lobes are well fit by a 59" Gaussian 
(dotted line). The near sidelobe level — 14 dB (—15 dB) of the 
type-1 (type-2) profile is consistent with optical crosstalk of 
2.5% (< 1%) between adjacent bolometers. 
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C. Band 

The expected bandpass for the system is calculated 
from the transmittance spectrum of the metal-mesh LP 
filters and the HFSS simulation of the waveguide, and 
cavity geometry. The spectral bandpass of sample chan- 
nels on type-1 and type-2 sub-arrays were measured using 
a Fourier transform spectrometer (FTS). The calculated 
and measured bandpass functions for both detector types 
are shown in Fig. 
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FIG. 15. Bandpass for a perfect cavity (dotted line) and 
that of type-1 (black) and type-2 (blue) cavity geome- 
tries. The bandpass of the perfect cavity is calculated 
from the mesh filter x waveguide transmittance. The type- 
1 and 2 detectors have curves for the FTS measured band- 
pass (solid lines with uncertainties) and the calculated mesh 
filter x waveguide x cavity transmittance (dashed lines). Esti- 
mated statistical measurement uncertainties are shown; these 
are not uniform for the measured type-2 spectrum due to 
post-measurement corrections for differences between the fil- 
ters in the test setup and those in the field. In addition, there 
are systematic uncertainties due to the optical coupling ge- 
ometry between the detectors and the FTS which may cause 
ripples in spectrum. The FTS measured results are normal- 
ized to the peak calculated transmittance of the correspond- 
ing detector type. The absorbed power of each configuration 
is proportional to the area under the curve. The measured 
type-1 bandpass is much narrower than expected, while the 
measured and calculated bandpass of the type-2 configuration 
are more similar. 



Each spectrum is characterized by an effective band- 
width, center frequency, and peak value. The effective 
bandwidth is given by 



A^ off = J F{v)du, 



(5) 



value of unity. The center frequency is 

/ F{v)v dv 
Vc = fF(v) dv' 

The expected loss in efficiency due to the type-1 (type- 
2) absorbing cavity is estimated by comparing the calcu- 
lated spectrum for the mesh filters, waveguide, and cav- 
ity to that assuming a perfect absorbing cavity (i.e, mesh 
filters and waveguide only) . We then define the cavity ef- 
ficiency, r] C av, as the ratio of the peak transmittances of 
these two calculated spectra. The cavity efficiency is 0.72 
and 0.78 for the type-1 and type-2 detectors, respectively. 







Calculated 


Measured 


Wafer geometry 


Tjcav 


Av cH (GHz) 


Av cH (GHz) 


Perfect cavity 


l 


37.6 




Type-1 


0.72 


32.3 


24.5 


Type-2 


0.78 


37.5 


33.5 



TABLE VI. Cavity efficiency and effective bandwidth calcu- 
lated for predicted and measured band pass spectra for the 
two detector types. 

The cavity efficiency as well as calculated and mea- 

Note 



sured effective bandwidths are given in Table VI 
that the FTS measurements are in relative units, so there 
is no direct measurement of Tj cav . The design bandwidth 
is 37.6 GHz. Taking into account the frequency depen- 
dence of simulated cavity absorption, the expected ef- 
fective bandwidth of the type-1 detectors is 32.3 GHz, 
while their measured Ai/ e ff is 24.5 GHz. The anoma- 
lously narrow band adversely affects the overall sensitiv- 



where F{y) is the bandpass function normalized to a peak 



ity of the type-1 detectors (Sec. VII H I. In contrast, the 
measured bandwidth of the type-2 detectors, 33.5 GHz, is 
much closer to the corresponding calculated bandwidth, 
37.5 GHz. While we do not fully understand the below- 
expected performance of the type-1 sub-arrays, the type- 
2 detector performance approaches that expected from 
simulations. 

Infrared leaks in the LP mesh filters can contribute 
significantly to the detector optical loading. We place 
upper limits on the high frequency out-of-band response 
of the instrument in the laboratory with measurements 
using a series of high-pass thick grille filters with cutoffs 
at 146 GHz, 206 GHz, and 305 GHz. The signal from a 
hot-chopped source sets an upper limit of 3x 10 -4 for the 
ratio of out-of-band to in-band response to an Rayleigh- 
Jeans (RJ) temperature source. 



D. Calibration 

The detector timestream data are recorded in counts 
from the ADC which are proportional to the bolometer 
current. We make a daily raster scan of an astronomical 
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source of known brightness to measure the conversion 
from ADC units to astronomical source flux. 

Mars is our most frequently used calibration source. 
The brightness temperature of Mars is calculated from 
the Rudy modeP^ scaled by a factor 1.052. The scaling 
factor is derived by comparing the Rudy model results 
at 93 GHz with Wilkinson Microwave Anisotropy Probe 
(WMAP) measurements of Mars at 93 GHz (additional 
details in Ref. 62). When Mars is not available, we ob- 
serve the RCW38 HII region, which we have calibrated 
at 150 GHz with Mars. 



1. Calibration Procedure 

We calculate a factor to convert ADC counts to 
changes in RJ temperature, ATrj, on each detector for 
each night of observation. The equivalent RJ tempera- 
ture change on the sky seen by detector n is 



AT R j n — 



AP 



S.n 



2k B Av cS ' 



(6) 



where ks is the Boltzmann constant, Aiy c g is given by 
Eq. (jjjj), and APg„ is the expected optical power dif- 
ference between the calibration source S and the CMB 
for a single-moded beam. The individual calibration 
scan maps from each detector are corrected for the fi- 
nite source size and fit to an elliptical Gaussian beam as 
described in Sec. |VIIB 1| Then, the calibration factor for 
each detector is 



AT, 



a T RJ ,n — 



RJ.n 



A,,, 



(7) 



where A n is the amplitude of the Gaussian fit to the 
uncalibrated source map. 

Determining the expected power from the source, 
APg yj, requires a measurement of the beam pattern. As 
detailed above (Sec. VII B 2), we cannot accurately mea- 



sure the beam solid angle for each beam. We, therefore, 
use the mean beam profile (Fig. 14 ) and assume a fixed 
ratio of the true beam solid angle to the best-fit Gaussian 
beam, 



n 



p r , dn 



p 



Here P n and G n are the normalized beam and elliptical 
best fit Gaussian power patterns, respectively, and p is 
the scaling factor between the two, for example, p — 
1.32 for type-1 detectors. The best-fit Gaussian beam 
solid angle for each channel is increased by this factor. 
With this approximation, the beam dilution factor for a 
compact source, such as Mars, is calculated as 



VS.n 



(8) 



where the integral in the numerator is taken over the 
solid angle of the calibration source. 



With the above form for r)s <n , APs jn is 

ao 2 f B(v,T s ) - B(v,T C mb) „, s , , n , 
APs,» = 7is, n c / 5 F{v)dv, (9) 



where B{v 7 T) is the Plank blackbody brightness spec- 
trum, Ts is the source temperature, and Tcmb is 2.73 K. 
For the beam measurement and calibration of August de- 
tailed here, Mars was used with Tg = 207 ± 3 K 

The calibration factor is calculated from equations 
([6]), Q, ([8]), and Finally, we apply an elevation 

dependent atmospheric opacity correction to each scan. 
The correction has the form 



X = e 



where r is the optical depth at zenith, e scan is the median 
elevation of the scan to be calibrated and e ca i is the ele- 
vation of the calibration scan. Section IVII D 21 describes 
the measurement of r. The correction factor is generally 
< 3%. 



2. Atmospheric Opacity 

Atmospheric emission is measured with skydips in 
which the telescope is scanned from zenith to low ele- 
vation on each night of observation. Modeling the atmo- 
sphere as a horizontal layer of emission produces a sky 
signal which will vary with observation elevation (e) as 



ATrj(e) = T atm (e 



-r csc(e) 



(10) 



where T atm is the temperature of the atmosphere and r 
is the optical depth at zenith. 

In fitting Eq. (10 1, r and T atm are largely degenerate. 



We lack reliable atmospheric temperature data, so we 
assume T atm — 273 K. Figure [16] shows five different 
skydips for a range of optical depths taken in August 
2007 along with model fits. In the fits, each channel is 
fit individually, but the plots show the median response 
at each elevation with a curve set by the median r. The 
skydips are well fit down to an elevation of 30° , indicating 
that the detector response remains linear over this range 
of opacities. 

In the August 2007 observing period, the optical depth 
at zenith ranged between 0.017 and 0.072, with a me- 
dian of 0.032. Using the APEX atmospheric transmit- 
tance calculator and the APEX-SZ spectral bandpass, 
we calculate a PWV of 1.4 mm for the median optical 
depth. We extrapolate the results to zero airmass and 
determine the background load contribution of the at- 
mosphere. There is an atmospheric load of ~ 1.3 pW on 
the detectors at 55° elevation when the optical depth is 
0.032. 



3. Calibration Stability 

The calibration uncertainty is determined by compar- 
ing the calibration from separate days. For seven succes- 
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FIG. 16. Skydip scans of varying opacity. We fit each scan 
according to Eq. ( 10 1 to find r, assuming an atmospheric 
temperature of 273 K. The curves shown here have, from top 
to bottom, t = 0.072, 0.046, 0.035, 0.025, 0.017. 



sive days, we fit the coadded Mars map to a Gaussian. 
Among the seven maps, there is a ±1.9% variation in 
fit amplitude after correcting for differences due to Mars 
temperature and atmospheric opacity. This variation is 
due to small changes in bolometer bias parameters from 
night to night. 

We typically focus the secondary mirror and perform 
only one large array calibration raster for each night of 
observation. Focus drifts change the relative calibration 
of detectors but have little effect on the average array 
calibration. Defocusing of beams in one region of the 
array is largely compensated by focusing the opposite 
region. 

The measured array responsivity can change during 
the night due to changes in the temperature of the detec- 
tor cold stage or atmospheric loading on the detectors. 
Short scans of calibration sources are used to monitor 
gain stability over the course of the night. In a sample 
series of Mars observations taken over 4 h, there was a 
±1.6% brightness variation among scan maps after cor- 
recting for elevation changes. Therefore, while the aver- 
age calibration on a given night is known to 1.9%, the 
average detector gain for an individual scan differs by an 
additional 1.6%. 



E. Pointing 

A pointing model was developed for APEX from op- 
tical observations of stars. The measured pointing error 
is ±2" and tracking accuracy is ±0.6" under stable at- 
mospheric conditions However, we frequently observe 
through sunset and sunrise during which temperature 
changes cause slow drifts in the pointing. Such drifts 
increase the effective beam size in maps coadded from 



multiple scans. During observations, we check the point- 
ing by scanning a planet or quasar every 1-2 h or when 
changing targets. The applied corrections have an rms 
variation of ~ 4", which we take to be the pointing un- 
certainty. 

Measurements of the positions of the individual beams 
relative to the center of the focal plane are very stable. 
The rms variation for Mars scans, which vary 20° in el- 
evation, is 1.1". No correction is necessary. The total 
pointing uncertainty for the beams is ±4.1". The correc- 
tion to the calibration from this uncertainty is less than 
0.5%, which is negligible. 



F. Instrument Efficiency 



The change in power measured in the bolometer, 
AP opt , and the expected power from a source, AP$ in 
Eq. ([9]), differ because of losses in the optical elements, 
loss of thermal energy in the spiderweb, imperfect cavity 
absorption, and uncertainty in the conversion of thermal 
to electrical power in the bolometer. We define an effi- 
ciency associated with each of these elements: r] optl r\boio-, 



and r\ T 



such that 



TJtot — Vopt ?]bolo ?]cav T]resp- 

The optical efficiency, rj opt , is calculated from the ex- 
pected reflection and scattering of optical elements. The 
transmittance of individual elements are listed in Ta- 



ble VII However, 34% of the throughput of each feed 
horn is truncated at the Lyot stop. Since this is an ele- 
ment of the optics design rather than a loss due an imper- 
fect optical element, we compensate for this loss in the ef- 
ficiency calculation. By this measure, a perfectly lossless 
optical path which includes a Lyot stop has rj op t = 1. The 
total transmittance of the optical elements other than the 
Lyot is r] op t = 0.62. 

The bolometer absorption efficiency, ry&o/oi is discussed 
in Sec. |IV| Calculation of the cavity efficiency, r] cav , is 
described in Sec. IvTTCl 

The responsivity efficiency, ry resp , can deviate from 
unity due to finite ETF loop gain as shown in Eq. Q. 
The bolometers are biased at a high temperature in the 
transition and are operated with a loop gain of 5-10. We 
estimate that the responsivity differs from the infinite 
loop gain limit by up to 15% but do not have an accu- 
rate measurement of the loop gain for all detectors. In 
addition, the measurement of V^ias deviates from the ac- 
tual bias voltage applied to the bolometer by up to 10% 
due to phase shifts in the cryostat wiring and residual 
impedance of the LC resonator. However, these two fac- 
tors have opposing effects on the power responsivity. As 
a rough approximation, we assume that r\ resp is unity. 

The cumulative efficiency with which the instrument 
converts changes in source power to electrical power in 
the detector, rjtot, is measured by the ratio of AP opt to 
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APs- Thus, the cumulative efficiency is given by 



Vtot 



Op l 



AP s {i-L Lvot y 



(ii) 



where L^ yot — 0.34. The median measured APEX-SZ 
efficiency is ritot = 0.31 for type-1 detectors and 0.36 
for type- 2. The difference in efficiency is predicted by 
simulations of the two backshort geometries. 



G. Optical Loading 

We measure the total background optical load on the 
bolometers by comparing the detector current as a func- 
tion of bias measured during observation with those of 
dark detectors. Dark voltage response curves were mea- 
sured in the laboratory before deployment. The median 
total optical load on the type-1 bolometers during the 
August 2007 observations is 6.1 pW. This corresponds to 
a loading temperature Trj = 44 K. 

The estimated contribution to the optical load for each 
optical element is listed in Table |VII| The atmospheric 
load is the measured value described in Sec. lVIID 21 The 
total estimated load is 6.5 pW, and estimated efficiency 
is 0.36 ( 0.24/(1 — L Lyot ) ). Both are greater than the 
measured values. The most likely source of error is r\ resp . 
An r] reS p of 0.9 would result in a calculated load and 
efficiency consistent with the measured values. 

Atmospheric loading is the largest single component 
of the background optical load. The sky and telescope 
contribute 40% of the total loading. Significant contribu- 
tions also come from the reimaging mirrors and the 60 K 
thermal filters. The three mirrors inside the cabin have 
high quality surfaces, but they scatter light inside the 
cabin to 300 K. Together, the mirrors contribute roughly 
35% of the total optical load. 



H. Noise and Sensitivity 

I. Instrument Noise 



Table VIII summarizes the expected median noise con- 



tributions from readout, detectors, and optical load for a 
single detector. Values used to calculate the noise terms 
are the median values of the array for type-1 (type-2) de- 
tectors: Vuas = 6.0 (3.3) iN, Puas = 22 (11) pW, and 
P opt = 6.1 (12) pW. 

The SQUID and warm readout electronics (readout 
amplifier, feedback resistor, and nulling resistor) produce 
broadband current noise which appears incoherently in 
both sidebands of the carrier. When the bolometer sig- 
nal is recovered, the sidebands are summed and the noise 
from both contribute. In addition, these signals do not 
pass through the LC filters so, when demodulated with 
the square- wave mixer, there are additional contributions 
from noise at odd harmonics of the carrier band. The to- 
tal noise is a factor n/2 larger than the nominal noise 



level of these components. We refer these noise cur- 
rent terms to a noise equivalent power (NEP) via the 
TES responsivity, Eq. Q. Note that the ratio of sig- 
nal to readout noise for this system, which uses alter- 
nating voltage bias and a square-wave demodulator, is 
within 10% of that for a constant voltage bias detector. 
This ratio would be unity if the demodulator in this sys- 
tem used a sine-wave mixer, as is the case for the digital 
multiplexer PS 

The bolometers have a noise contribution from the 
fluctuation of thermal carriers in the bolometer thermal 
link. The thermal carrier noise has the form given in 
Table [VTTT] The factor 7 accounts for a temperature gra- 
dient along the linkP^l For a normal metal link at the 
detector operating temperatures, 7 = 0.58. The Johnson 
noise from the TES and bias resistor also appear incoher- 
ently in both sidebands of the carrier. These signals are 
filtered by the LC filters so there are no additional har- 
monic contributions. The total noise is a factor \[2 larger 
than the nominal Johnson noise of a resistive element. 
When operated at high loop gains, ETF suppresses noise 
current though the bolometer, including Johnson noise. 



Table VIII includes the Johnson noise at the full value, 
which should be taken as a worst case. 

We include two terms for the photon noise from back- 
ground loading: shot noise and "bunching" noised The 
bunching term takes into account the boson nature of 
photons and has an uncertain magnitude which we pa- 
rameterize as £ € [0,1].^ The expected total detector 
white noise is in the range 64-75 and 65-92 aW/vTlz for 
type-1 and type-2 detectors, respectively. 

In the type-1 detectors, the total readout noise, 
bolometer thermal noise, and photon noise are all compa- 
rable. The narrow bandwidth and 0.72 cavity efficiency 
of the type-1 detector reduce the incident optical power. 
As a result, the thermal conductance, chosen based on 
the expected load, is unnecessarily high. The high G re- 
sults in a high thermal carrier noise and increases the 
Vbias required to maintain the detectors in the supercon- 
ducting transition. Thus, all the non-optical noise terms 
are elevated. The electrical bias power required to hold 
the bolometers in the superconducting transition is 3.5 
times the absorbed optical power. 

Type-2 detectors move toward background limited 
performance with increased bandwidth and optimized 
bolometer properties. The wider bandwidth results in 
larger P opt and increased photon noise. The lower G 
decreases the required Vbias and P e i ec which decreases 



all non-photon noise contributions (see Table III). The 
lower R n decreases the required V i as which results in 
lower readout noise. In the type-2 detectors, applied bias 
power and absorbed optical power are roughly equal. 



2. Detector Sensitivity 

The achieved bolometer NEP is the white noise level 
in the signal band after atmospheric noise removal 
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TABLE VII. Optical loading and loss for individual optical elements. T c is the physical temperature of each element, e is 
the emissivity, L a is the spillover/scattering loss of each element, T s is the temperature of the spillover/scattered radiation 
absorber, r> e is the estimated transmittance/efficiency of each element, and P opt is the contribution to total background optical 
power absorbed by the bolometer. Results are tabulated for type-1 bolometers. The total load corresponds to an RJ loading 
temperature of 44 K. 



(Sec. VIII). The signal band is the white noise re- 



gion between 1/f noise and the bolometer response roll 
off. After atmospheric noise removal, the bolometer 
timestreams have 1 // noise with a knee at 1-2 Hz which 
is dominated by the detector readout. We take 3 Hz to 
be the minimum frequency which is safely above the 1//. 
The bolometer time constant (Tabic III ) filters signal fre- 
quencies / > l/(2/TT opt ). The resultant signal band is 
3-14 Hz and 3-8 Hz for the type-1 and type-2 detectors, 
respectively. When observing, the telescope scan speed 
is selected so the modulated sky signal appears in the 
signal band (Sec. |VIII A ). 



Tabic |IX| lists the median sensitivity of the bolome- 
ters in flux and temperature units, while Fig. [TT] shows 
a histogram of the pixel sensitivities of the type-1 array 
and a type-2 sub-array. The values are calculated from 
bolometer channels which remain after analysis pipeline 
data cuts. 

The type-1 detectors have a median NEP of 
97 aW/VHz, median NET of 890 /^K C mb%/s and me- 
dian NEy of 3.5xl0~ 4 *Js. Here, y is the dimensionless 
Comptonization parameter, proportional to the cluster 



optical depth. The type-1 NEP is well above the maxi- 
mum expected noise of 75 aW/VHz. Approximately one 
third of the channels have an NEP in the expected range. 
These are the channels with lower bias frequencies. As 
described in Sec.|VJ there are frequency dependent phase 
shifts between the carrier and output demodulator in the 
readout system. The phase shifts result in the sensitiv- 
ity being dependent on bias frequency which is not ac- 
counted for in Table IVHll and which elevates the detector 
noise above the expected level for high bias frequencies. 
The lower optical load on the type-1 detectors enhances 
the frequency dependent excess noise because the non- 
optical noise terms are more significant. This frequency 
dependent effect increases the variance in detector sensi- 
tivity of Fig. [17| The long tail of high NET detectors is 
due to a subset of detectors with low optical efficiency. 

The NEP of the type-2 detectors, 87 aW/VTlz, is 
within the expected range. The type-2 detectors also 
show elevated noise with increasing bias frequency, but 
the effect is smaller, and the median NEP falls within 
the rage of expected values. The type-2 detectors have 
a median NET of 530 /iKcmbv^ an d median NEy of 
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Noise Source Equation Type-1 NEP (aW/VHz) Type-2 NEP (aW/^Hz) 

Readout 

Warm electronics 3.7 X 7r/2 [pA/-\/Hz] X Vm */V2 25 14 

SQUID 2.5 x tt/2 [pA/VSi] x Vuas/V2 17 9.2 

Bias resistor Johnson noise 1.4 x \/2[pA/Vnz] x Vbi as /y/2 8.5 7.9 

Bo/ometers 



Thermal carrier noise \J^k B T£G 36 33 

Johnson noise V4fc s T c P Mas 24 19 

Photon shot noise xf2hvP opt 35 50 

Photon bunching noise yJ^P^ pt / Av e s 39 66 

Total (£ = 0) 64 65 

Total (£ = 1) 75 92 



TABLE VIII. Summary of contributions to NEP. Median values of the array are used to calculate the expected noise contri- 
butions from readout, detectors and optical loading. The photon bunching noise is listed as an upper limit corresponding to 



Sensitivity 
NEP (aW/VHz) 
NEFD (mJy^i) 
NET GuKrjVs) 
NET (^KcMBx/i) 
NEy (xl0~ 4 y/a) 



Type-1 
97 
42 
500 
890 
3.5 



Type-2 
87 
24 
300 
530 
2.2 



TABLE IX. Observed median noise and sensitivity per chan- 
nel. NEP is the measured noise equivalent power in the detec- 
tor. NEFD is the noise equivalent flux density, measuring the 
sky-signal sensitivity. The NET values are the noise equiv- 
alent temperature referred to a source at the RJ limit and 
the CMB temperature. NEy is the noise equivalent y, the 
dimensionless Comptonization parameter. 



2.2 xl0~ 4 s/s. Overall, the type-2 sub-arrays have me- 
dian detector sensitivities a factor 1.7 better than the 
type-1 design. This improvement is consistent with the 
wider measured bandwidth and higher cavity efficiency 
of these detectors. 



VIII. OBSERVATIONS AND ANALYSIS 

We operate the receiver on a 24 h cycle, in line with 
telescope sun avoidance periods and the sorption fridge 
hold time. Each night, we check the secondary mirror 
focus position, perform 1-2 skydips to characterize the 
atmosphere (Sec. VII D 2), and make periodic pointing 
checks (Sec. VII E). Approximately 10% of our observa- 
tion time is used for flux calibration, including at least 
one full raster of a known source for calibration and beam 
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FIG. 17. Histogram of individual channel sensitivities in the 
array. Sensitivities are measured in the signal band of each 
detector type after analysis pipeline atmospheric removal and 
data cuts. NET values are plotted for the 135 type-1 channels 
(dark gray) which remain after all data cuts during an obser- 
vation of the Bullet cluster. Also included are 20 channels 
from a type-2 sub-array (blue) measured in a separate obser- 
vation. The type-2 histogram does not obscure any type-1 
values. 



measurement and periodic scans of secondary calibrators 
to monitor gain stability through the night (Sec. VII D 3 1 . 
With these tasks, target switching, and telescope down 
time, roughly 12 h each night are spent on science tar- 
gets. 
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A. Scan Strategy 



B. Cluster Map Making 



A key goal of APEX-SZ is the development of a scaling 
relationship between cluster mass and SZE flux. To this 
end, we have focused on targeted observations of known 
clusters, both relaxed and evolving systems, over a range 
of redshifts. 

We use circular drift scans to efficiently map a single 
cluster target within the wide APEX-SZ FOV. In the cir- 
cular drift scan, the telescope makes circular scans which 
are stationary in azimuth and elevation. The source 
drifts across the field with the Earth's rotation. After the 
source drifts across the array, the telescope is reset for the 
new source position. Figure [18] illustrates a typical scan 
of twenty 12' diameter circles with a period of 5 s. This 
circular scan mode offers a few advantages. The scan has 
a continuous, low acceleration, therefore, no data are dis- 
carded due to high acceleration turnarounds. While low 
acceleration, the scan speed is fast enough (~ 8 arcsec/s) 
that the signal from a typical cluster appears well above 
the 1/f knee of readout and atmospheric noise (Fig.[l9|. 
In addition, the scan is fixed in AZ/EL coordinates allow- 
ing us to search for ground signal or other contaminants 
by summing data from multiple circles. 
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After completing each scan, the data are recorded as 
280 readout channel timestreams sampled at 100 Hz 
with telescope pointing data interpolated at the same 
timesteps. The data reduction pipeline used for generat- 
ing cluster maps from raw timestreams generally follows 
a simple sequence: cuts to remove low quality data and 
optically dead channels, filtering to remove atmospheric 
noise, cuts to remove data which have residual noise af- 
ter filtering, and binning filtered data into maps. In ad- 
dition to channel cuts (Sec. VII A), timestream data are 



rejected due to step glitches from electrical interference, 
high telescope acceleration at the start and end of scans, 
and noisy segments in otherwise good channels. After 
all cuts, roughly 40% of the recorded 280 channel data 
contributes to the map. 

Atmospheric noise mitigation is based on a low-order 
time-domain polynomial, which removes slow drifts in 
the timestreams, and low-order spatial polynomial, which 
removes signals correlated among bolometer channels ac- 
cording to their relative position in the array. Figure [l9| 
shows noise spectra for a single bolometer during a scan 
of the Bullet cluster before and after atmospheric noise 
removal. After filtering the large scale correlated noise 
from the detector time ordered data, the 1/f knee of the 
remaining noise is lowered to ~1 Hz. Since we can recover 
the predicted white noise level in the low bias frequency 
channels, we are confident that we have successfully re- 
moved the atmospheric noise in the signal band. 
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FIG. 18. The circular drift scan pattern used for mapping 
individual clusters. The telescope scans in 12' diameter circles 
with a 5 s period (black line). The circle center is stationary 
in AZ/EL, but tracks in RA/DEC. After performing twenty 
circles, the telescope is reset for a subsequent circular drift 
scan. Also shown are the FOV of the bolometer array at the 
initial pointing (blue dashed line) and the source position (x). 



FIG. 19. Measured power spectral density (PSD) for a type- 
1 bolometer channel during a single scan made in weather 
typical of the August 2007 observing run. The dashed line 
shows the noise spectrum as measured. Its 1/f knee is 10 Hz. 
The solid line shows the same data after atmospheric noise 
removal by the analysis pipeline. The 1/f knee of the data 
after atmospheric removal is at 1.4 Hz, below the signal band. 



After atmospheric noise removal, the timestream of 
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each channel is binned in a map with the channel AZ /EL 
offset and flux calibration applied. These are coadded to- 
gether with inverse variance weighting to make a map for 
each scan. Finally, the individual scan maps are coadded 
with inverse variance weighting to make the final map of 
the cluster. 

Reference \§21 details the data reduction and map- 
making procedures used to create the SZE map of the 
Bullet cluster shown in Fig. [20] The Bullet cluster (IE 
0657-56) is a large cluster with total mass ~ 10 15 M & and 
gas temperature ~ 1.4 x 10 s K. It consists of two merg- 
ing sub-clusters: the eponymous bullet passing through a 
larger cluster. APEX-SZ observations of the Bullet were 
made in August 2007. Because of the large size of the 
Bullet, the data for this map are processed with the clus- 
ter masked in order to preserve signal at the expense of 
increased noise in the map. After data processing and 
final smoothing with a 1' Gaussian kernal, the resulting 
map has 85" resolution with 55 y,K rms noise, resulting 
in a 23fi detection in the central 1' region of the map. 



The map of the Bullet cluster is made from the reduc- 
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FIG. 20. SZE map of the Bullet cluster based on 6.5 hr of 
observation. Contours are spaces at 100 /iK. The map has 
85" resolution with 55 /j.K rm3 noise. 



tion of 234 circular drift scans like that of Fig. 18 for a 
total of 6.5 h of integration time. The integration time 
is typical of the APEX-SZ cluster maps, which require 
of order 10 h of integration per cluster. Because of the 
large FOV of APEX-SZ, a significant fraction of the data 
is taken well away from the region of interest around the 
cluster. The equivalent of 3.6 h is spent integrating the 
central 20' x 20' area of the map. 

The analysis and model fitting of the Bullet cluster 
in Ref. [62] is performed on lower noise map produced 
without masking the cluster during noise processing. It 
has 27 iiK rms noise over the central 20' x 20', for a 
mapping speed of 1.0xl0~ 4 deg 2 /VKc MB /hr or, equiva- 
lently, 4.7xl0 -2 deg 2 /mJy 2 /hr. There are a few factors 
which contribute to the discrepancy between the mea- 
sured and calculated mapping speed. Residual correlated 
noise likely remains in the map because the extended size 
of the Bullet requires relatively gentle noise filtering to 
avoid removing cluster signal. Additionally, there are 
several high noise detectors which pass data cuts. These 
are included in the integration time, but, when inverse 
variance weighted, make little contribution to the final 
map. 

Extrapolating from the Bullet results, 30 h of obser- 
vation on a single target would result in a map with 
13 /iKqmb noise in the central area. This integration 
time is a typical deep integration, and constitutes a 
rough practical limit to map sensitivity achievable with 
APEX-SZ. The sensitivity limit is compatible with our 
goal of calibrating large SZE surveys, which will map 
> 1000 deg 2 to a depth of 10-20 /jK CM b- 

The XLSSU J022145.2-034614 cluster is the lowest 
mass cluster detected by APEX-SZ, with an x-ray tem- 
perature of 4.8 ±0.6 keV and estimated mass A/500 — ^ x 
1O 14 /i _1 M P A 5ct detection in a map with 12 /iK CM B 
noise per V pixel, the cluster represents the effective clus- 
ter mass sensitivity of the instrument. 



IX. CONCLUSIONS 



C. Mapping Speed 



by 



The expected mapping speed of the instrument is given 



5' 



•^V bolo^beam^ldata 

NET 2 



(12) 



Nboio is the number of detectors, £lbeam is the beam 
solid angle and rjdata is fraction of collected data which 
contributes to the map. Applying the measured de- 
tector properties, with Nboio — 280, f2 



beam 



=1 amin", 

and rjobs = 0.4, we use equation (|12[) to calculate a 



nominal mapping speed for the type-1 detector array of 
1.4xl0- 4 deg 2 /fjK 2 CMB /hx. 



Wc have presented the design and described the perfor- 
mance of APEX-SZ, an instrument primarily designed to 
observe the SZE in galaxy clusters at 150 GHz. APEX- 
SZ is the result of several years of development and em- 
ploys robust, scalable technologies in a modular design 
which have already been adapted and significantly im- 
proved for use in other instruments. The focal plane is 
a lithographed TES bolometer array with a 280-channel 
frequency multiplexed SQUID readout system. The ar- 
ray is cooled by a pulse-tube cooler and helium sorption 
refrigerator to an operating temperature of 280 mK. 

Data from observations were used to analyze detec- 
tor yield, beams, bandwidth, optical efficiency, optical 
loading, and sensitivity. The detector beams show good 
uniformity across the array and a low level of crosstalk 
between adjacent channels. The atmosphere and tele- 
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scope are the primary optical load. Despite a high peak 
optical efficiency, the sensitivity of the type-1 detectors 
is limited by a narrow bandwidth. The median sen- 
sitivity per type-1 channel is 890 [iKcmb\/s (NEy of 
3.5xl0 -4 i/i). In 2009, we implemented a type-2 de- 
tector sub-array with an improved cavity design with a 
A/4 backshort and a median sensitivity of 530 (jKcmbV^ 
(NEy of 2.2xl0~ 4 y/s). 

APEX-SZ has observed for a total of 875 h since 2007. 
We have produced high signal to noise images of x-ray 
selected clusters and produced low noise maps of fields 
chosen to overlap with external data sets. The cluster 
map of Fig. [20] and modeling of the Bullet cluster are 
presented in detail in Ref. 1621 A power spectrum analysis 
of our observations of the XMM-LS^3 field is given m 
Ref. Ell In Ref. E3 APEX-SZ and x-ray observations of 
Abell 2163 are used to deproject the thermal structure of 
the intra-cluster medium. Reference [70] details a similar 
analysis of Abell 2204. 

Observations were primarily focused on a sample of x- 
ray selected clusters spanning a wide range of mass and 
dynamical state. APEX-SZ has been used to image 48 
clusters. These observations will be used to study the 
scaling of the SZE signal with masses derived from weak 
lensing, x-rays, or optical richness for a sample of clusters. 
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